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Abstract 

We study the critical velocity of a Bose-condensed gas in a moving one-dimensional (ID) optical 
lattice potential at finite temperatures. Solving the Gross-Pitaeavskii equation and the Bogoliubov 
equations, within the Popov approximation, we calculate the Bogoliubov excitations with varying 
lattice velocity. From the condition of the negative excitation energy, we determine the critical 
velocity as a function of the lattice depth and the temperature. We find that the critical velocity 
decreases rapidly with increasing the temperature; this result is consistent with the experimental 
observations. Moreover, the critical velocity shows a rapid decrease with increasing lattice depth. 
This tendency is much more significant than in the previous works ignoring the effect of thermal 
excitations in the radial direction. 

PACS numbers: 03.75.Kk, 03.75.Lm, 67.25.de 
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I. INTRODUCTION 



The occurrence of energetic and dynamical instabilities in a Bose-Einstein condensate 
(BEC) moving in a periodic (optical lattice) potential is an interesting problem from the 
conceptual viewpoint, since it involves basic properties of superfluidity. Experimentally, en- 
ergetic instabilities of BECs moving through optical lattices have been observed both in the 



weak lattice and tight-binding regimes ]. According to Landau's argument 



4J, the 



energetic instability of the superfluid state is attributed to the appearance of negative exci- 
tation energy (which is known as Landau instability). Once a superfluid velocity is beyond 
its unique critical velocity, the current suffers friction leading to the decay of superfluidity. 
We note that the actual occurrence of energetic instability requires some dissipative mecha- 
nism, which leads the system to a lower energy state. Some theoretical studies showed that 
the thermal component plays a crucial role, receiving the energy emitted by the condensate 
during the breakdown process . Experimentally, energetic instability is only observed 
at finite temperatures, consistent with the theory of Refs 5, f|. Thus, in order to provide 
a quantitative account for the critical velocity relevant to the experimental data, it is im- 
portant to consider thermal component at finite temperatures. However, most theoretical 
studies on the Landau instability of a Bose gas in an optical lattice used the zero-temperature 
GP equation [?], Isj]. Moreover, they concentrated on the first Bloch band using the Bose- 
Hubbard model, and have ignored the effect of radial excitations. It will turn out that one 
should include thermal excitation in the radial direction to obtain quantitative results. 

In this paper, we study the critical velocity of current carrying condensate in a ID optical 
lattice at finite temperatures. For this purpose, we calculate the Bogoliubov excitation 
energy in the framework of Hatrree-Fock-Bogoliubov-popov (HFB-Popov) approximation in 
a periodic potential, with explicitly including the effect of the radial excitations. In Sec. II, 
we derive a quasi-lD model of the Gross-Pitaeavskii equation and the Bogoliubov equations 
for the Bose gas moving in a ID optical lattice. Using the HFB-Popov approximation, 
we solve these equations to calculate the Bogoliubov excitation spectrum. In Sec. Ill, we 
calculate the critical velocity of current carrying condensate in a optical lattice, which will 
be determined from the condition of the negative excitation energy. We obtain the critical 
velocity as a function of the lattice depth with a fixed temperature, and compare it with the 
experimental data The magnitude of the critical velocity is found to be in reasonable 
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agreement of the experimental data [lj. We will also show that the critical velocity decreases 
with increasing lattice depth, which is consistent with the experimental result 1|. We also 
calculate the temperature dependence of the critical velocity with a fixed lattice depth. The 
critical velocity will be found to drop very rapidly with increasing temperature. 



II. QUASI ID MODELING OF A CURRENT CARRYING CONDENSATE 



We consider a Bose condensed gas in a combined potential of highly elongated harmonic 
trap and ID optical lattice. Our system is described by the following Hamiltonian : 



H 



h 2 9 



+|«t(r)i.t(r)tf.(r)tf.(r) 



ext I 



V>(r) 



(1) 



where g 



4wh 2 a 



is the coupling constant determined by the s-wave scattering length a. The 



external potential V ext is given by V ext (r) = Vt rap (r) + V op (z). More explicitly, we consider 
a cylindrical condensate, which is radially confined by the harmonic potential V tra , p (r) = 
Y [ui^x 2 + y 2 )], and is subject to the periodic potential V op (z) = sEr cos 2 (^z). Here s is a 
dimensionless parameter describing the intensity of the laser beam creating the ID lattice 
in units of the recoil energy En = 4p-, where ^ = ^ is fixed by the wavelength A of the 
laser beam. In this paper, we neglect the harmonic trap potential along the z-direction. we 
consider a highly anisotropic cigar-shaped harmonic trap potential. In order to take into 
account this quasi-lD situation, we expand the field operator in terms of the radial wave 
function 



(2) 



where (p a (x,y) is the eigenfunction of the radial part of the single-particle Hamiltonian 

h 2 



___ vi+ ^ +!/2) 



1pa(x,y) = £ai>a(x,y), 



(3) 



which satisfy the orthonormality condition f dxdyip^(x, y)ijjp{x } y) = 5 a p. Here a = (n x , n y ) 
is the index of the single-particle state with the eigenvalue e^,^) = hjj^{n x + n y + 1), where 
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u± is the trap frequency is the radial direction. Inserting Eq. (|2j) into Eq. ([I]) and using 
Eq. (J3j) , we obtain 

h 2 d 2 



H = J dztjj a (z) 



2m dz 2 



Vop(z) + € c 



aa'/3/3' 



dzipli'l'lpp^a'- 



(4) 



The renormalized coupling constant is defined by g aa 'pp' = 9 J dxdy<l>a<f>p<l>p , <f>a'- Following 
the procedure described in Ref. 11], we separate out the condensate wavefunction from 
the field operator as ^ a = (ip a ) + ip a = & a + i/) a , where $ a = (ip a ) is the condensate 
wavefunction and ip a is the noncondensate field operator. The Popov approximation neglects 
the anomalous correlation 11]. Within the HFB-Popov approximation, we obtain the 

generalized Gross-Pitaevskii (GP) equation 

h 2 d 2 



io| 



2m dz 2 + + 



®a+Yl 9aa'pp' (®*p® p'$a> + 2$? ^J^ a '}) • (5) 



a'PP' 



Solving the GP equation (5) at T = (setting (^'fyj = ) using the trap frequencies 
relevant to the experiment l| (see the paragraph below Eq. (jl]) ), we find that fi <C Hlj± and 
|$ a | 2 /|$o| 2 ^ 10~ 6 {ol 7^ 0), where we have denoted the lowest radial state as a = = (0, 0). 
Thus the contribution from higher radial modes to the condensate wave function is negligible 
small. For this reason, we will henceforth approximate $ a ~ $<5 aj o- Taking the usual 



Bogoliubov transformations for the noncondensate, ip a (z,t) = £V fc 
obtain the coupled Bogoliubov equations, 



u jka a jk - v* ka a] k 



we 



L a v,j a + y] Q i 



L a Vj a + X^a' 



(2s£' n o + g aa >pp>nppi) u ja > - g%n v ja > 



(2#£ra + gaa'pp'nppi) v ja > - g%n u ja > 



~E jUjoti 



(6) 



where we have introduce the operator L a = —^-§^2+ K P (^) + £ a — fJ> 



We have also 

introduced the simplified notations n Q (z) = \&(z)\ 2 , g", = g aa '0Q,npp' = yftp'&P'J- As noted 
above, we only include the lowest mode (a = 0) in the condensation wavefunction $. Sums 
over the repeated indices [3,(3' are implied in Eq. (jSJ). These equations define the quasi- 
particle excitation energies Ej and the quasi-particle amplitudes Uj a and Vj a . Using the 
solutions of Eq. ([6]), one can obtain the noncondensate density from h = ^ a v, where 



n 



aP 



{u ja Ujp + VjaVjp) N(Ej) + V ja V jl3 



with N(Ej) = 1/ [exp(pEj) - 1]. 
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We expand the condensate wavefunction in terms of the reciprocal lattice vector G, 
^( z ) = Z~2n e^ g+nG ^ z C n . Here we consider the situation where the condensate is moving with 
the velocity hq. We also expand the Bogoliubov quasiparticle amplitudes u ka and v ka in 



terms of the reciprocal lattice vector, u ka = u kn e 



a p i(nG+k)z 



Vkc 



J2 n v% n e i( - nG+k > } where 



k is the Bloch wavevector (or quasimomentum) of excitations. Then, GP equation and 
Bogoliubov equations become 

Vfi 2 



U (9+nG) ' 



+ eo — \i H — — ^ j C n H — — ^ (C n+ i + C n -\) + g C* m ,C m C n 



+m'—m 



a'/3 jkW 



(7) 
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V jkn-1 



jkn+l'-l 



vf kl vf kl , (Nf k + 1) 



iifi' 



jkW 



V 



jkn+l'-l 



/3 mm' 







^ m' m 




jkn+m'+m 




(8) 



jkn— m' — m 

The dispersion relation for the condensate moving through the optical lattice is obtained 
by solving the Bogoliubov equations in Eq. (jHJ). Solving the coupled equations ([2j) and (jHJ), 
we self-consistently determine the excitations spectrum Ej and the condensate fraction at 
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finite temperatures. Our calculation procedure closely follows that summarized in Refs. 
10] . Throughout this paper we use the following parameters of the experiment of Ref. []J: 
m ( 87 Rb)=1.44xl0" 25 kg, uj z /2tt = 9.0 Hz, u ± /2ix = 92 Hz, scattering length a = 5.82 nm 
and the wavelength of the optical lattice A=795 nm. We fixed the number of atoms per 
lattice site as J_^, 2 ndz = 200. 

In Fig. [T]we show an excitation spectrum for s = 1, T = 20nK (~ 0.14T C ) and q = 0. 
The q = case corresponds to the condensate at rest. We only show the first two Bloch 
bands (j =1 and 2). For each band, we plot the first 20 radial branches a = 0, 1,2, 19. 
This spectrum is consistent with Modugno's results for T = [l^ . 



In Fig. [21 we plot the condensate fraction N c /N as a function of the temperature for 
various values of the lattice depth s. In the actual calculation, we must introduce cutoff 
in the summation over the radial modes(a) and the reciprocal lattice vector (n). We chose 
sufficiently large cutoff so that the results are cutoff-independent. It is crucial here that 
the cutoff is strongly temperature-dependent, because the thermal occupation of elementary 
excitations strongly depends on temperature. 

As shown in Fig. [21 the transition temperature deceases with increasing lattice depth. In 
this paper, we focus on this weak lattice region s < 5. In a stronger lattice s > 10, one can 
use the tight binding approximation [l3]. In contrast, a full HFB-Popov formalism presented 
in this paper must be used for a shallow lattice potential, say, s < 5. In this paper, we only 
present the numerical results for s < 5, but we expect the results to smoothly continue for 
s > 5. 
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FIG. 1: Excitation spectrum of a cylindrical condensate at rest in the optical lattice (quasimomen- 
tum of the condensate q=0) as a function of the quasimomentum k of the excitations, for lattice 
depth s = 1. The first two Bloch bands are shown and for each band we plot the first 20 radial 
branches. The lowest branch in each band corresponds to axial excitations with no radial nodes. 



III. THE CRITICAL VELOCITY OF CURRENT CARRYING CONDENSATE 

We now consider the excitation spectra of the lowest phononic branch (a = 0) for s — 1 
at T=20nK(~ 0.14T C ) for different values of the condensate velocity (Fig. [3]). By increasing 
q, the dispersion law becomes asymmetric; the slope of the phononic branch decreases for 
the mode propagating opposite to the condensate current (k < 0). By further increasing q, 
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T[nK) 

FIG. 2: The condensate fraction as a function of temperature for different values of the lattice 
depth s = 1, 3, 5 
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FIG. 3: The excitation spectrum of the current carrying condensate with a fixed lattice depth 
s=l and temperature T = 20nK(~ 0.14T C ), for different values of the condensate momentum q 
(q/G=0, 0.01, 0.02). We only plot the lowest radial branch (corresponding to axial excitations 
with no radial nodes). 

the slope of the spectrum at k — > —0 becomes zero at a certain critical value q c . Beyond 
q c , low energy excitations become negative. Therefore, superfluidity breaks down at q c . For 
comparison, we also calculate the Bogoliubov sound velocity c for q = for different values 
of the temperatures T. We obtained sound velocity c from the slope of the Bogoliubov 
excitation spectrum at k — > — 0. In Fig. HJ we plot the temperature dependence of the 
critical velocity q c and the sound velocity c with fixing the lattice depth as s = 1. It is 
clear that critical velocity is always lower than the sound velocity. We also find that critical 
velocity decreases rapidly with increasing temperature. The reduction of the critical velocity 
is more significant than that of the sound velocity. The reduction of both these velocities 
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are due to the reduction of the condensate fraction with increasing temperature, which is 
larger for the condensate with a finite q than the condensate at rest. This tendency is much 
more significant than in the previous works ignoring the effect of thermal excitations in the 
radial direction. 

We now compare our result with the experiment of Ref. In the experiment with 35% 
thermal fraction, a strong reduction of the number of atoms was observed even with an 
extremely low velocity, which indicates that critical velocity at finite T is much lower than 
the sound velocity. Ref. [l| attributed this low critical velocity to the harmonic confinement 
in the z direction. In fact, Ref. [jjj] pointed out that the lowering of both the critical velocity 
and sound velocity is caused by inhomogeneous density profile of the harmonically trapped 
BEC. Another possibility of the observed extremely low critical velocity is a sudden onset of 
dissipation in the low density tail of the condensate cloud Q, [jjj]. On the other hand, from 
our calculation, the temperature corresponding to 35% thermal fraction is estimated as T ~ 
70.5nk. As shown in Fig.|4]at this temperature, the critical velocity has a strong temperature 
dependence, and the critical velocity at 70nK is much lower than the sound velocity, i.e. 
hq c /c(T = 0) < 1. Therefore, even ignoring the trap confinement in the z direction, we 
find that the critical velocity becomes extremely low at the temperature relevant to the 
experiment In any event, in order to make quantitative comparison between our results 
and experiments, it will be necessary to remove the effect of inhomogeneous density of 
the harmonically trapped BEC (for example, using the technique used in the experiment 
of Ref. [14I ). We note that it is crucial to include the thermal excitations in the radial 
direction in order to obtain the strong reduction of the critical velocity. 

In Fig. [5j we plot the Bogoliubov sound velocity c and the critical velocity q c as a function 
of lattice depth s with fixing the temperature as T = 20nK. We find that the critical velocity 
rapidly decreases with increasing lattice depth s. Fig. [5] also shows that the critical velocity 
is always lower than the sound velocity c. The reductions of both the critical and sound 
velocities are due to the increase of the effective mass m* with increasing lattice depth. As 
shown in Refs. 0, [si, Q, , m* depends on the condensate momentum, and is larger for 
finite q. This gives rise to the significant change in the critical velocity. 
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FIG. 4: Temperature dependence of Bo- 
goliubov sound velocity (q = 0) and the 
critical velocity with a fixed lattice depth 

8=1. 

IV. CONCLUSION 
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FIG. 5: Dependence of Bogoliubov 
sound velocity and the critical velocity 
on the lattice depth s with a fixed tem- 
perature T = 20nK. 



In this paper, we studied the temperature dependence of the critical velocity of current 
carrying condensate in a ID optical lattice, with explicitly including the effect of the radial 
excitations. In the trap geometry we considered, the condensate wavefunction can be treated 
only with the lowest radial mode, but the radial excitations for thermal cloud cannot be 
neglected. 

Within the Popov approximation, we calculated the temperature dependence of the Bo- 
goliubov excitations with varying lattice velocity. From the condition of the negative exci- 
tation energy, we determined the critical velocity as a function of the lattice depth and the 
temperature. For comparison, we also calculated the sound velocity of condensate at rest 
in a ID optical lattice potentials. As shown in Fig. HI the critical velocity rapidly decrease 
with increasing temperature, much faster than the sound velocity. The difference between 
the sound velocity and critical velocity was observed in the experiment of Ref . [lj . Although 
the direct comparison of our results with the experimental data is difficult because of in- 
homogeneous density profile in the z direction, the critical velocity is shown to be strongly 
reduced at the temperature relevant to the experiment. 
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